(t  j 


NAYORD  REPORT 


*)■ 

OX 

m 

C3^ 

01 

i 

a:) 

«ca.. 

UNDERWATER  EXPLOSION  PHENOMENA: 
THE  PARAMETERS  OF  MIGRATING  BUBBLES 


t  js: 

Uo 


^  cn 

■' ' 


m 

mHaW 


12  OCTOBER  1962 


H£liSAS£D  to  JtSTXA 

HAVAL  OHDNAUCE  LABOKATOKt 
^^ithout  restrictions 

□  for  . elease  to  Military  and  Government 
Asenci  s  Only. 

□  Approval  by  BuWeps  required  for  releiiae 
to  contraotors. 

□  Approval  by  BuWeps  required  tor  all 
subsequent  release* 


U.  S.  NAVAL  ORDNANCE  LABORATOR' 

^  WHITE  OAK.  MARYLANi: 


;  C  \  '• 

’  -'rif— " 


Best 

Available 

Copy 


FOR  ERRATA 


THE  FOLLOWING  PAGES  ARE  CHANGE! 

« 

TO  BASIC  DOCUMENT 


FOR  ERRATA 


A  9^  y:(y 


THE  FOLLOWING  PAGES  ARE  CHANGES 
« 

TO  BASIC  DOCUMENT 


mORD  I;l85 


' 


'  y 


ft  ' 
y 


a: 


< 


UNDERWATER  EXPLOSION  PHENOMENA:  THE  PARAMETERS 
OF  MIO^TING  BUBBLES 


Ly 

Haas  G, 


ABSTRACT:  The  migration  of  underwater  explosion  hubbies  caused  by 
(*  loyancy)  affects  the  energy  of  the  pulsation  as  well  as  period  er 
radius  in  the  second  aad  subsequent  cycles*  Experimental  data  on  p 
explosions  in  various  depths  are  analyzed^  and  the  bubble  energy  is 
for  live  cycles  of  the  oscillation  as  a  fianction  of  the  strength  oi 
The  result  is  given  in  dimensionless  form  which  perinits  the  calcula 
t  periods,  maximum  radii,  and  migration  for  a  wide  range  of  conditior 

An  energy  balance  shows  the  surprising  resul*::  that  the  bubble  energ 
•  second  cycle  increases  with  increasing  intensity  of  migration  imtil 

is  reached  at  a  condition  of  strong  mlgiation*  Beyond  this  point, 
decreases  again*  This  gain  is  found  to  be  due  to  the  decrease  of  t 
radiated  by  the  bubble  pulse* 
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NAVORD  REPORT  4lB5  12 

Uin)E2KWATER  EXPLOSION  PHEN(»ffir.:  THE  PA  A^2n•1^2^S  OF  MIGRATING  BUBB 

In  this  r^ort  the  attempt  is  made  tc  further  the  understanding  o 
hydrodynamic  processes  associated  vith  migrating  e:q)loslon  bubble 
results  are  obtained  from  an  evaluation  of  experimental  obse:*vati 
the  use  of  involved  theoretical  calculations.  The  imports 
calculations  is  s  graph  showing  the  bubble  energy  ratios  foj  ^ 
cycles  of  the  oscilM^cion,  Combined  with  simple  equations  this  i 
permits  the  calculation  of  the  bubble  parameters  for  almost  all  p 
ixnpjrtant  conditions  of  common  eaqploslves  detonated  under  water, 
ccaning  NOLTR  62-1B4  utilizes  the  results  of  this  paper  and  presen 
method  for  expedient  calculations. 

Preliminary  results  of  this  study  have  been  reported  in  previous 
Althou^  this  final  version  shows  lit-ule  differences,  the  previou 
the  bubble  energy  ratios  must  be  considered  to  be  siiperseded. 

This  report  is  part  of  a  comprehensive  study  on  the  behavior  of  e 
bubbles,  work  has  been  caiT-*ed  out  under  Task  No,  301-664/43 
and  RUME-3-S-000/212-1/WFC08-10-004  PA  002. 
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HAVORD  REPORT  4lfi5  12  October  I962 

UNDERWATER  EXPLOSION  PHEl«»iEHA:  THE  PARAMETERS  OF  MICaSATING  K  iBLRS 

'i  ■^^is  report  the  atterapii  is  made  to  Tirther  the  understanding  of  the  various 
.  rodyn&mic  processes  associated  vlth  migrating  explosion  bubbles*  Revealing 
results  are  obtained  from  an  evaluation  of  experimental  observations  vithout 
the  use  of  involved  theoretical  calculations*  The  isiportant  outcome  of  the 
calculations  i.*?  a  graph  showing  the  bubble  energy  ratios  for  the  first  five 
cycl  i  of  the  oscillation*  Combired  with  sJiple  equations  this  information 
permits  the  calculation  of  the  bubole  parameters  for  almost  all  practically 
important  conditions  of  coimon  explosives  detonated  under  water*  The  forth- 
ccnuxig  NOLTR  62-184  utilizes  the  results  of  t^is  paper  and  presents  a  simple 
method  for  e3q)6dient  calculations* 

rellminary  results  of  this  study  have  been  reported  in  previous  papers* 
/hou^  th  3  final  version  shows  little  differences,  the  previous  data  on 
*.he  bubble  ener  *^  ratios  must  be  considered  to  be  superseded* 

This  report  is  part  of  a  compreheiisive  study  on  the  behavior  of  explosion 
bubbles*  'Ihe  work  i  xs  been  carr-*ed  out  under  Task  No.  301-664/43003/01040 
and  RUME-3-E-000/212-1/WFC08-10-004  PA  002* 
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!•  DmiOrUCTION 

There  is  an  extensive  literatisre  on  the  p'xlsatioa  o  huVules  produced  by 
underwater  explosions,  listed  in  references  (a)  and  (t  Almost  all  of  this 
vork  refers  to  non«migrating  bubbles,  i«e«  it  is  aasu^  ed  that  *;he  center  of 
the  bubble  remains  stationary*  Tjils  assumption  is  valid  for  small  charges 
e^iploded  in  great  depth,  e*g*  1  Ib  in  500  ft*  The  bubbles  of  leirg'^r  explosions 
move  upward  because  of  their  buoyancy  and  thereby  change  the  characteristics 
of  the  pulsation  narameters*  Little  is  known  today  about  the  fundamental 
behavicr  of  such  oubbles  in  the  sectvrd  and  the  following  cycles* 

A  theoretical  treatment  of  this  problem  is  difficult  for  two  reasons: 

(a)  Migrating  bubbles  change  their  shape*  When  contracting,  the  Lower  inter¬ 
face  of  the  originally  spherical  bubble  moves  faster  inward  than  the  ipper 
interface*  At  an  intermediate  moment,  the  cross  section  of  the  bubble 
resembles  that  of  a  kidney*  Later,  the  upper  and  lover  interfaces  collide 
and  the  bubble  beccmes  a  torus*  Upon  re-expansion,  the  spherical  shape  is 
rou^ly  restored,  but  energy  has  been  dissipated  by  the  inpact  of  the  inter¬ 
faces*  (b)  Near  the  bubble  minimum,  the  gas-water  interface  becomes  unstable* 

It  tends  to  dissolve  into  a  water  spray  \diich  is  projected  into  the  interior 
of  the  bubble*  This  brings  forth  a  cooling  of  the  esqplosion  gases  and,  thus, 
again  a  dissipation  of  energy* 

Both  of  these  dissipative  processes  are  difficialt  to  account  for  theoret¬ 
ically*  A  further  portion  of  the  bubble  energy  is  radiated  by  the  bubble 
pulse*  These  energies  are  not  available  for  the  subsequent  cycles  of  the 
bubble  pul  ‘tion*  A  knowledge  of  the  r gaining  energy  is  of  prime  importance 
in  any  quantitative  calculation  of  the  later  bubble  phenomena. 

A  rather  crude  approach  is  used  in  this  paper  to  find  the  bubble  energy 
in  each  cycle  cf  the  oscillation*  The  analysis  is  based  on  the  change  of 
the  bubble  period  caiised  by  the  various  degrees  of  migration*  Strong  migration 
carries  the  bubble  into  a  shallower  depth*  As  a  consequence  the  period  of  the 
next  oscillation  is  Increased  whea  cooipared  with  that  of  a  non-migratixig  bubble* 
To  carry  out  the  analysis,  the  magnitude  of  the  b^^bble  migration  xoust  be  known* 
This  process  has  been  e:q>erimeatally  studied  by  means  of  acoustic  ranging  of 


*  The  list  of  references  is  found  at  the  end  of  this  report. 
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the  biibble  miuixmnn*  Further  evidence  is  obtained  in  this  paper  frcxn  the 
observation  the  time  at  'which  the  bubble  breaks  the  water  surface*  Once 
a  biibble  migrates  into  the  proximity  of  the  water  surface,  the  emission  of 
the  bubble  pulse  ceases*  The  evidence  'vAiether  or  not  a  bubble  pulse  is 
observed  before  the  breakthrou^  yields  information  on  the  d'  «ch  of  the 
preceding  bubble  maximum* 
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n.  EQUATION  FOR  THE  BUBBI£  ENERGY  RATIO 


The  period  of  the  bubble  oscillation  in  the  n-tli  cy-  e  is  given  by 


0.374  t  (1-0  . 

2n 


The  symbols  are  esqplalned  In  the  list  at  the  end  of  this  report#  The  first 
portion  of  (l)  lb  given  in  referenco  (b),  where  rQ  is  denoted  by  e#  The  last 
term  in  the  parenthesis  accounts  for  the  effect  of  the  vater  surface  uhich 
shortens  the  period  This  term  vlU  be  discussed  in  paragraph  HI#  The  effect 
of  the  bottom  of  the  sea  is  not  included  in  (l)>  because  the  tests  i«hlch  ve 
Mill  analyze  were  made  in  vater  sufficiently  deep  as  to  mate  this  effect 
negligibly  small# 

Forming  the  ratio  of  the  n-th  to  the  first  period,  ve  have 

(2)  vr.  rT7^- 

In  this  equation,  the  dimensionless  period  t  has  been  cancelled  vhlch  amounts 
to  the  tacit  assimption  that  it  is  the  same  for  all  cycles  of  the  pulsation# 

On  the  basis  of  the  classic  bubble  theory,  this  is  valid  if  the  adiabatic 
exponent  ganma  of  the  explosion  gases  is  between  1#2  and  1#3  or  so#  Then,  t 
is  essentir'^ly  constant  and  independent  of  the  aa^litude  of  the  pulsation,  as 
seen  in  Figure  2  of  reference  (b)#  The  gamna  of  most  of  the  common  e^losives 
is  within  this  range#  For  ,ther  values  of  gamma,  the  above  assumption  is  only 
approximately  valid# 

The  maxisium  bubble  radius  in  the  n-th  cycle  is,  reference  (b): 

.  r  QW  .1/3 

(3)  V  ■  (-!_) 

n 

A  surface  correction  is  not  necessary  for  the  maximum  radius,  reference  (a) 

Chapter  8#  Assuming  that  is  the  same  for  all  cycles,  ve  can  eliminate 

A_  in  (2)  by  means  of 
nn 
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The  assim^tlon  of  an  equal  is  net  as  good  an  apprex  jatla^  as  tliat 
of  eq.ual  t*  Since  the  significance  of  the  aaxiinum  radii  in  the  later  pilses 
is  less  critical  than  that  of  the  periods,  the  acc\rracy  is  acceptable  for 
our  analysis* 

Re^^arrangesnent  of  (2)  fields  the  fol3  jig  quadratic  for  the  cube  root 
of  the  bubble  energy  •‘•itio 


(5) 


Qf 


0. 


In  order  to  find  eaqpreasions  for  the  ratios  of  the  hyirootr.tic  head 
and  the  depth  nigration  of  the  bubble  must  be  Imovn*  We  refex' 

to  the  center  of  the  n-th  bubble  maximum  and  ve  set  for  the  rise  betveen 
tyro  such  points 


(6) 


n 


IThis  fonmila  is  ccmmonly  used  for  the  migration  betveen  the  point  of  esq^losion 
and  the  first  minimum*  It  represents  an  approximation  of  Taylor’s  fo*”x"ula 
derived  in  reference  (c)* 

If  the  assumption  is  made  that  the  bubble  remains  sphericed,  theoretical 
calculations  yield  for  the  coefficient  C*  the  value  132*  (For  the  rise 
betveen  two  successive  maxima,  C*  in  (6)  vould  be  about  tvice  as  large*) 
Conipari8<»i  of  the  theoretical  value  vith  that  obtained  frm  measured  bubble 
migrationo  shovB  the  fonper  about  h(A  too  hi^di*  Figure  1  illustrates  this 
discrepancy  and  also  illustrates  the  great  scatter  of  the  data* 

JL  ^2  2 

The  magnitude  w  '‘’/Z  in  (6)  is  relaw^d  to  the  Froude  Nuniber  F  vhich 
for  our  purposes  is  best  defined  by 
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(7) 


52;  “-C  —  y  • 


ThllS 


(8) 


^  V 
Z  -  ^ 


-3/2 


For  the  we 
Into  the  form 


in  conjunction  vlth  the  quadratic  (3};  (6)  can  be  brought 


Then 


(10) 


“-1  D-D 

I  ^ 

>1 


and 


(n)  J 


n-l 

I 

J-1 


iLliaia 


The  computation  of  the  reduced  migration  in  the  n»th  cycle  requires 
only  the  knovledge  of  provides  aH  relations 

necessary  for  the  evaliSticm  of  (3;*  The  energy  ratios  r^/r^  can  be  computed 
if  the  foil  ring  infomatloa  is  given: 

(a)  numerical  values  foi  0  and  a 

T  A%i  D 

(b)  ^  as  a  function  of  ‘=—  and  ~  • 

1  1  1 

The  next  paragraphs  deal  vlth  these  data* 
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ni.  THE  SURFACE  OORRECTION  TERM 


The  surface  correction  term  in  (l)  is  given  in  the  literature  in  various 
forms.  Reference  (d)  lists  the  same  form  ob  (l),  but  footnote  12,  peigc  3^2, 
of  reference  (a)  shows  the  surface  correction  as 


(12) 


1  - 


(Reference  (•)  enroneously  shows  a  ^/6  power  for  Z, )  For  the  case  of  a  free 
water  surface  and  infinitely  deep  water,  F'(x)  is  unity*  With  the  introduc¬ 
tion  of  (3)  and  a^^  0,92  (which  is  appropriate  for  our  conditions),  one 

obtains  or  *  0*21^, 

Essentially  the  sajae  value  is  obtained  for  the  coefficient  occurring  In 
the  siirface  correction  term  for  the  migration  (6)  for  which  we  have  used  the 
same  symbol  a*  From  equation  (8«67)  of  reference  (a),  this  magnitude  is  found 
to  be  0,2*  (There  are  two  misprints  in  reference  (a)  at  this  place:  In 
equation  (8*67)  the  factor  0*2  is  omitted.  In  the  preceding  eqmtion,  the 
factor  should  read  0,4  and  not  0,2),  Reference  (d)  quotes  a  =«  l/5  for  the 
siirface  correction  of  the  migration. 

These  values  as  well  as  the  form  of  the  surface  correction  tern  are  first 
order  approximations.  Figure  2  ^ows  the  period  constants  versus  depth  of 
es^iosion  corrected  with  the  uee  of  various  coefficients  a.  These  data  are 
from  a  test  series  which  contributed  the  most  ii^portant  valiies  to  our  evalua¬ 
tion,  It  is  seen  that  cr  *  0,1  makes  the  period  constant  independent  of  depth, 
i,e,  eliminates  the  surface  effect,  whereas  a  «  0,2  ’"overcorrects".  Actually, 
cr  is  not  a  constant,  but  depends  on  A^/O  as  well  as  figure  8,21  of 

reference  (a)  shows  fair  agreement  for  300  Ib  TUT  charges  and  Figure  8,20  good 
agreement  fc'  0,66  lb  tetryl  charges,  if  the  periods  are  corrected  with  or  «  0,2, 
But,  a  closer  inspection  again  reveals  an  overcorrection,  if  shallow  firing 
coriitions  are  excluded.  It  s*  ims  that  a  *  0,1  is  appropriate  for  such  cases 
where  the  bubble  is  not  too  close  to  the  surface,  as  in  Figure  2,  Where 
<  0,5,  For  larger  values  of  for  shallower  e^losioos,  a  value  of 

a  larger  than  0,1  is  niedec  or  ccooplete  correction.  This  shows  that  the 
sinple  form  of  the  surface  correction  term  in  (l)  is  not  sufficient  for  a 
precise  description  of  this  effect.  For  the  purpose  cf  uils  paper  the  ccanpH- 
cations  of  a  more  elaborate  relationship  would  not  be  worthwhile.  In  the 
majority  of  cases  which  are  of  interest  here,  the  bubble  is  so  far  away  from 
the  water  surface,  that  the  value  or  »  ,1  is  appli-'^able.  For  the  Instances 
discussed  In  paragraph  IV,  where  the  bul  *16  break-throu^  is  considered,  it 
was  found  that  ot  had  little  effect  on  the  l(‘;c*ation  of  the  bwbble  and 

on  the  times  of  the  break-throu^.  Thus,  it  seems  seife  to  use  cr  *  0,1  for  the 
purposes  of  this  paper. 
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IV.  OHE  MIGRATION  COEFFICIENT 


In  Figure  1  es^erimental  data  on  the  bubble  migration  up  to  the  first 
bubble  Mnifflum  are  coii5>iled.  The  reduced  migration  AZ/F-  is  >l'jtted  versus 
100  The  data  stem  from  various  experimental  seri'  >  a  variety 

of  charge  weights  and  explosives.  Most  of  the  experjjner  si  points  were  taken 
from  a  compilation  in  reference  (f).  Onlj*^  typical  poin\  5  are  shctwn  in  Figure  1 
in  order  to  avoid  an  overloading  of  the  graph.  AH  values  are  reduced  to  an 
explosive  hrving  the  properties  of  TNT  with  the  use  of  the  appropriate  equiva- 
lent  factors.  Furthermore,  all  data  are  corrected  to  free  field  conditions, 
i.e.  an  infinite  medium.  The  correction  was  made  as  described  in  Paragraph  III 
of  this  paper  using  the  factor  a  =  0.3*  The  results  of  the  290  lb  TNT  tests 
are  not  very  certain,  because  they  were  carried  out;;  in  water  oi  limited  d^th. 
The  method  used  to  coi'^ect  for  the  effect  of  the  bottom  consisted  of  including 
2 

the  term  ~  cif  into  the  parentheses  of  (6),  \Aiere  H  is  the  distance  be¬ 

tween  bottom  and  point  of  explosion.  This  correction  is  probably  not  very 
accurate,  but  the  only  one  known  today. 

The  migration  was  measured  by  the  sound  ranging  meunod.  The  point  from 
which  the  bubble  pulse  is  emitted  was  located  by  means  of  triangulation  from 
a  vertical  string  carrying  several  pressure  gages.  One  of  the  difficulties 
of  this  method  arises  from  the  fact  that  the  pulse  is  not  emitted  from  a  mathe¬ 
matical  point,  but  from  the  bubble  surface  which  at  the  moment  of  the  minimum 
still  has  a  considerable  size.  The  great  scatter  of  data  shown  in  this  figure 
Illustrates  the  diffic\iltieB  of  this  measurement  and,  therefore,  the  rather 
approximate  nature  of  the  information. 

Figure  1  also  shows  the  result  of  numerical  calculations  using  Taylor’s 
migration  fc  -nula,  reference  (c).  The  points  shown  have  been  obtained  by  rather 
laborious  computations  made  during  World  War  II  by  varioxis  British  agencies 
quoted  in  reference  (f).  Re.  erence  (d)  stated  that  on  this  basis  the  rise  is 
proportional  to  Kennard,  reference  (e),has  noticed  that  the  line 

drawn  through  these  theoretical  points  can  be  equally  well  represented  by  a 
simple  relationship  ^ich  in  our  notation  corresponds  to 


(13) 


The  experimentally  observed  migration  .  r-.n  then  be  represented  by 


(14) 


c* 


y 
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vhere  C*  is  between  80  and  90  depending  on  the  weight  one  desires  to  gire  the 
migration  at  shallow  or  greater  depths  of  e^^losion  respectively.  Actually 
Figure  1  suggests  a  sli^tly  different  functional  dependency  of  W  and  Z  from 
that  of  the  above  eqtuatlon.  It  is  dubious  whether  or  not  the  '  .perimentcQ. 
evidence  is  stiff icientV  accurate  to  establish  such  a  functl'^'^.,vl  relationship# 
In  view  of  the  scatter  and  the  difficulties  of  meetsurlng  mi  ration,  the  form 
ilk)  appears  to  be  adeqmte#  Also,  it  seems  that  either  of  ihe  above  quoted 
values  of  the  C*  fits  the  data  equally  well  and  can  be  used  with  equal  con¬ 
fidence. 

For  the  purpose  the  present  analysis  the  rn5.gration  between  successive 
bubble  is  needed  and  not  that  between  the  point  of  explosion  and  the 

bubble  minimum  ^ich  is  shown  in  Figure  1.  Ihe  theory  of  reference  (g) 

(which  deals  with  sphericax  bubbles)  gives  the  following  picture  about  these 
migration  terms: 

Migration  between 

(a)  point  of  esg^losion  and  first  maximum 

2 

(15)  [in  2  -  1/2  ]  , 

(b)  point  of  explosion  and  first  minimum 

g  T  ^ 

(16)  Dq  -  In  ^  c^, 

(c)  first  minlnum  and  second  maytmiin 

n.  2 

(17)  [  In  2  -  1/2  +  In  4  J 

(d)  first  laeodmum  and  second  naximum 

(18)  D^-  Dg  -  I  In  4  c^  +  Tg®  In  4  Cg 

+  (t/  -  )(ln  2  -  1/2-)] 

The  magnitude  c  depends  on  the  ratio  of  maxlimna  to  TP.Inlmum  bubble  radius: 

(19)  c  ^ 
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As  before,  the  subscripts  1  and  2  refer  to  the  first  and  second  cycle  of  the 
bubble  pulsation  respectively* 

On  the  basis  of  this  theory,  the  migration  between  the  tvo  maxima  is  Just 
twice  the  migration  between  the  point  of  explosion  and  t^ie  .irst  minimum,  if 
the  period  and  the  coefficient  c  are  equal  in  bolJh  cycle**.  /vltb»>ugji  this 
result  gives  a  vsiluable  hint,  the  approximations  made  1  .  (iC),  (Ll),  and  (6) 
are  apparent:  (a)  Migration  betwesn  the  point  of  explosion  and  f-he  first 
bubble  maximum  is  neglected  in  (lO)  and  (ll)  for  simplicity*  Relation  (l^) 
which  is  probably  rather  accurate,  since  all  bubbles  remain  spherical  during 
the  first  esg^anslon,  indicates  that  this  migration  term  can  indeed  be  neglected* 
It  is  small  in  ccr^srison  with  the  other  terms  if  the  firing  conditions  are 
such  tha'«  sevexal  cycles  of  bubble  puiration  occur*  It  is  a  poor  approxima¬ 
tion  for  large  charges  e^gploded  so  shallow  that  the  breakthrou^  at  the  water 
surface  occurs  after  ihe  first  minimum*  But,  such  conditions  are  not  the 
subject  of  our  study*  (b)  Since  neither  T  nor  c  are  consistently  equal  for 
successive  cycles,  the  factor  2  is  not  necessarily  applicable  in  (6)*  But,  nei¬ 
ther  is  any  constant  factor,  since  the  migration  between  tvo  maxima  depends  on 
the  parameters  of  both  cycles*  It  is  the  siinplicity  of  the  enalysis  as  well 
8LS  the  present  lack  of  any  better  information  which  Justifies  approximation 
(6)*  In  view  of  these  discussions,  it  appears  to  be  desirable  to  obtain  an 
Independent  check  for  the  migration  formula  (6)  and  the  factor  tt 

Such  an  additional  evidence  of  the  bubble  Mgration  can  be  obtained 
using  the  analysis  developed  in  this  paper  and  considering  the  time  at  which 
the  migrating  bubble  breaks  through  the  water  surface*  This  event  can  be 
conveniently  recorded  by  means  of  photography  of  the  surface  phenomena*  A 
pertinent  result  is  listed  in  Table  I* 

TABLE  I 


Charge  Weight;  1590  lb  TIB?  Equivalent 


Firing  Depth 

130 

140 

150  ft 

Time  of  Biibble  Break- 

through  Observed 
Kineratographicaily 

2.50 

3.24 

3*25  sec 

Time  of  Latest 

Bubble  Pulse 

2.55 

3.32 

2*90  sec 

Observed 

3rd  pulse 

4th  pulse 

4th  pulse 

Time  counts  from  moment  of  detonation* 
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The  sli^t  discrepancies  hetveen  the  observed  bubble  pulses  and  the 
tliaes  of  thi*  bubble  breakthrough  are  probably  experimental  errors*  One 
should  expect  the  time  of  breakthrou^  to  be  somevhat  larger  than  the  time 
of  the  pulse*  For  our  purpose  the  magnitude  of  the  time  intervals  observed 
is  not  ijai^>ortant,  but  only  the  correlation  of  these  events* 

Figures  3a  to  3c  show  the  position  and  size  of  the  bubbi  calcu* 

lated  with  the  coefficients  C  «  3*7,  3*5,  and  3*2  for  the  coLiltions  of 
Table  I*  Bbqplosion  bubbles  are  spherical  to  the  first  maxinium  only*  The 
circles  referring  to  the  subsequent  cycles  in  Figures  3a  to  c  must  be  luider** 
stood  as  Ideailized  measures  of  size*  Such  bubble  shapes  are  not  well  defined 
and  can  be  considered  spherical  in  a  crude  way  only*  The  nujzibers  shown  near 
the  bubble  x'efer  to  e  cycle  of  the  pulsation* 

At  130  ft  firing  depth,  the  center  of  the  fourth  bubble  TnAifinnim  occurs 
according  to  these  calcuijatlons  either  above  or  so  close  to  the  water  surface 
that  a  fourth  bubble  pulse  could  not  have,  been  emitted*  This  is  in  agreement 
with  the  experimental  evidence  which  shows  that  the  breakthrou^  must  havq 
occurred  tnmediately  after  the  3rd  bubble  minimum*  It  must  be  visualized 
that  the  bubble  center  imps  at  this  moment  rapidly  froa  the  position  "3"  to 
the  position  ”4",  Figure  3a  shows  that  for  C  3*7  and  3*5  a  considerable 
part  of  the  fourth  bubble  maxlmtim  is  above  the  water  surface*  This  would 
result  in  the  observed  surface  disturbances*  On  the  other  hand  the  bubble 
is  too  deep  at  the  third  cycle  for  C  >  3*2  in  ordtr  to  produce  such  surface 
effects*  This  shows  that  C  must  be  larger  than  3*2'.* 

At  140  ft  firing  depth,  a  fourth  bubble  pulse  was  observed*  For  C  *  3«7 
the  bubble  is  too  shallow  In  the  fourth  cycle  in  order  to  produce  a  pulse* 

It  is  generadly  assumed  that  mce  the  depth  of  the  bubble  center  is  less  than 
90^  of  the  maxiimim  radius  no  bubble  pulse  is  emitted*  (This  depth  is  called 
the  venting  depth*  We  prefer  the  term  **blow-ln"*)  Although  this  evidence 
is  establishe'*  for  tas  first  bubble  maximum  only,  it  is,  at  least  approxi¬ 
mately,  applicable  to  the  later  cycles  also*  It  ttims  out  that  in  this  case 
C  »  3*5  is  the  largest  v8Q.ue  fo  lAiich  a  fourth  bubble  pulse  can  be  expected* 

We  have  here  a  rather  sensitive  criterion  for  the  migration  coefficient  C* 

Figure  3c  shows  the  case  of  150  ft  firing  depth*  The  calculations  vtsing 

any  of  the  three  values  for  C  are  cce^tible  with  the  ejqperlmental  evidence* 

One  may  argue  that  C  «  3*2  again  appears  low  in  view  of  the  evidence  that  the 
breakthrough  essentially  coincided  with  the  fourth  bubble  ainimuip* 

It  is  obvious  that  the  data  available  are  not  sufficient  for  an  umubiguous 
determination  of  the  migration  parameta'  .  However,  it  is  significant  that  the 
value  C  «  3.*>  idilch  is  the  largest  value  conristent  wixh  the  evidence  of  130  ft 
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firing  depth  and  iftilch  holds  for  the  third  cycle  Is  conqpatible  vlth  the 
acoustically  loeasured  migrations  shown  In  Flgiire  1  ^Ich  hold  for  the  first 
cycle.  C*  and  C  are  interrelated  by 

(13)  C*  -  C. 

On  this  basis  C*  »  8o  corresponds  to  C  «  l«7d5*  s  refers  bo  migrations 
up  to  the  minimum.  The  value  needed  here  would  be  about  twice  as  large^  hence 
3*57* 


For  the  practical  calculations  it  was  decided  to  use  the  rounded  value 

C  -  3o5. 

It  is  realized  tha;  this  value  cannot  claim  a  high  degree  of  acc\nracy.  It  is 
consistent  with  the  experimental  evidence  available  today^  but  because  of  the 
great  scatter  of  the  data  a  consldexable  uncertainty  remains. 

In  the  following  analysis,  the  value  of  C  is  someifcat  less  critical  than 
it  mi^t  be  escpected.  The  bubble  energies  in  the  various  cycxes  of  the 
pulsation  are,  of  course,  dependent  on  C«  But,  when  these  bubble  energies  are 
later  used  to  calculate  bubble  parameters,  some  of  the  uncertainties  c(xmected 
with  C  will  be  eliminated.  For  instance,  the  periods  calculated  by  this  method 
will  €Llmost  exactly  reflect  the  input  periods.  !Qie  migration  of  the  bubble 
and  the  position  of  the  bubble  center  of  the  various  bubble  maxima  aie  more 
sensitive  to  the  migration  parameter.  But,  these  represent  the  best  informa** 
tion  available  today. 

A  sensitive  method  to  check  the  validity  of  our  method  is  the  observa¬ 
tion  of  the  migrating  bubble  in  a  hi^  gravity  tank.  A  preliminary  study 
resulted  in  excellent  atgreement.  A  comprehensive  test  series  aimed  at  a 
thorou^  check  of  the  migration  parajoeter  is  planned. 
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V,  EXEEEBEmL  INPUT 

In  addition  to  the  siuc^ace  correction  factor  or  and  the  bubble  migration 
coefficient  0,  information  Is  needed  on  the  period  ratios  'Hie  form  of 

equation  (^)  does  not  require  the  Imowledge  of  the  periods  by  ^delves.  Al8o» 
the  specific  firing  conditions  axe  not  needed  for  the  evaluatif''  or  ohls  equa¬ 
tion.  It  is  sufficient  to  3mow  the  ratio  ^  r«itlo 

The  nature  of  this  input  makes  It  possible  to  utilize  test  results  frooi  e3g;>lo- 
slve  ch8u:ges  having  different  e3g;>lo8lve  material,  different  charge  weights,  and 
different  firing  conditions.  Ihe  bulk  of  the  Information  stems  from  a  test 
series  listed  In  reference  (h).  (l>enoted  as  M-serles  in  Figures  4  and  ^•) 

Other  data  (C-serles)  are  from  the  files  of  the  E  Department  of  the  Naval  Ord- 
neuice  Laborasory. 

Figures  4  and  5  show  the  period  ratios  and  plotted  versus 

Aj^/Z^.  !I5ie  corresponding  values  of  given  in'^Figure  4  as  curves. 

The  limiting  values  for  great  depth  ->  O)  are  ^2^^  "  “ 

0.565. 


Information  on  the  fourth  bubble  period  is  sparse  and  uncertain.  The 
curve  in  Fi{>cire  6  Is  based  on  the  following  'uhree  vedues: 


TABLE  II 

v^ 

D  /Z 

1  1 

T, /r 

4  1 

0 

1.0 

0,51 

0.15 

0.81 

1.17 

0.166 

0.82 

1.34 

The  period  ratios  of  non-xaigratlng  bubbles  (A^^/Z^  ^  O)  show  a  slight 
variatloi  for  different  e:<plosive8,  reference  (j).”^e  v'llues  chosen  refer  to 
esqplosives  which  mosz  closely  resemble  those  employed  in  the  M-  and  C-^erles: 
minol  for  the  second  and  torpex  for  the  third  cycle.  The  ratio  for  xhe  fourth 
cycle  Is  estimated  from  the  TKf  result,  reference  (l),  on  the  basis  of  the 
trend  e^diibited  by  the  different  e3q>lodlves  in  the  preceding  cycles. 

A  err  >  estimate  of  r^/r]j^  can  be  beised  on  the  period  ratio  of  non-migrating 
TNT-bubbles,  T  /Tj^  ■  0.927,  reference  (i\  In  absence  of  any  tetter  information 
the  constant  v^ue  r^/r^^  «  0,80  ml^t  he  iised  as  a  rou^  approximation. 
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VI,  THE  EHERG^  OF  MI(2l/lTINa  BUBBLES 

The  result  of  the  ceJ-culations  outlined  above  is  shov^'  isi  Figure  6,  -where 
the  energy  fractions  plotted  versus  the  reduc  a  miration.  The 

energy  r  QW  refers  to  the  bubble  energy  at  the  instan*^  ^.f  xii^  i-th  bubble 
maxiinum*  It  Is  veil  known  that  the  major  portion  of  i  lis  energy  is  potential 
energy  stored  in  the  water  during  the  biibble  expansion#  %e  remainder  is  the 
internal  energy  of  the  gas#  Kinetic  energy  and  the  energy  of  migration  are 
negligibly  small  at  this  mesoent# 

At  each  b^^b‘L  ^  i  minimum  a  reduction  of  the  bubble  energy  takes  place# 

The  enerrv  ^tch  is  lost  at  this  poinw  essentially  con5)ri8es  two  terms: 

(a)  eneity  acoustically  radiated  by  the  bubble  pulse  and  (b)  dissipated 
energy#  Both  of  these  energy  terns  depend  on  the  degree  of  bubble  migration# 
Strong  migration  reduces  the  aunplitude  of  the  bubble  pulse^  thus  reduces  the 
energy  loss  due  to  the  acoustic  radiation# 

For  non-migrating  bubbles  the  energy  dissipation  -t  the  bubble  sLlnlmum 
is  probably  a  consequence  of  the  Taylor  instability  of  the  bubble  interface# 
With  Increasizig  migration  and  with  the  consequent  decrease  of  the  excess 
pressiure  In  the  bubble  near  its  znlcimum,  the  instability  decreaises#  Thus^ 
the  intensity  of  the  spray  projected  into  txie  bubble  interior  is  ailso  de« 
creased  and  so  is  the  energy  dissipation#  However,  migration  causes  the 
inversion  of  the  bubble  and  am  impinging  of  the  xQper  and  lower  btabble  inter¬ 
faces#  The  watei*  haxoner  amd  the  spray  formation  connected  with  it  are 
probably  the  alternative  mechanisns  of  energy  dissipation#  These  take  over 
with  increasing  Intensity  as  the  effects  of  the  Taylor  instability  decrease# 

These  considerations  iUustrate  the  important  role  of  migration  in  the 
partition  of  these  energy  t^Tms#  In  fact,  our  analysis  clearly  shows  that 
the  bubble  energy  fraction  depends  on  the  strength  of  the  migration#  It  also 
suggests  the  term  which  Is  used  in  Figure  6  as  abb^.issa# 

Figure  6  shows  the  energy  beOance  for  the  second  cycle  at  the  left  hand 
side#  The  bubble  energy  ratio  plotted  as  obtained  ftrom  our  analysis 

and  the  input  discussed  in  Paragropn  V#  This  ratio  refers  to  the  energy  of 
the  bubble  at  the  second  maximum#  Since  this  energy  is  supplied  frem  the 
bubble  energy  of  the  fix*st  cycle,  t^e  term  1  -  r^/r  r^resents  the  energy 
decrease  yiblch  occurs  at  the  first  ^/Ubble  mixilmum#  *^6  ^ergy  of  cUe  pulse 
%ftii;:h  is  emitted  at  the  end  of  the  first  and  the  beginning  of  the  second  cycle 
is  shewn  ais  a  dashed  curve#  (This  curve  is  a  rather  crude  estimate  made  for 
the  sake  of  illustration#)  The  two  cxirves  divide  the  plot  into  three  bwds 
the  hei^tsof  xAxich  represent  (a)  the  bubble  energy,  (b)  dissipated  energy,  and 
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(c)  the  pulse  energy*  In  the  dimensionless  diagram  shovn^  these  three  energies 
add  up  to  unity  which  is  to  say  that  the  hubhle  energy  of  the  first  cycle  spHts 
into  these  '*’^ree  terms* 

An  interesting  result  is  that  the  bubble  energy  of  the  second  :yole 
initially  increases  with  increasing  migration,  if  the  other  peurajs  .^ers  of 
the  ea^losion  are  held  constant*  In  view  of  the  preceding  coot  this 

is  to  be  ejected*  It  is  also  seen  that  the  dissipated  energ;  is  roughly 
constant  over  a  considerable  range  of  xrlgratlon  intensity*  Zlus  is  surprisix^, 
since  two  different  mechanisms  of  energy  dissipation  are  probably  involved, 
both  of  Which  depend  on  migration*  Althou^  one  of  them  decreases  and  the 
other  increases  with  increasing  migration,  a  roughly"  constant  sum  of  these 
factors,  as  suggested  by  Figure  6,  was  not  necessarily  to  be  expected* 

For  very  strong  migration,  the  bubble  inversion  and  the  consequent 
energy  dissipation  appear  to  be  so  violent  that  a  decrease  of  the  bubble 
energy  in  the  second  cycle  results*  For  such  conditions  the  energy  radiated 
by  the  bubble  pulse  is  practically  negligible* 

The  most  is^Kurtant  outcome  of  these  calculations  is  that  the  energy  ratios 
given  in  Figure  6  permit  the  calculation  of  the  parameters  of  niigrat^ng  btlbbles 
for  almost  all  explosion  conditions  of  practical  interest*  (Nuclear  explosions 
are,  of  course,  not  included*)  calculations  can  be  extended  with  reasonable 
confidence  up  to  the  end  of  the  fourth  cycle*  Estimates  for  the  fifth  cycle  are 
possible  with  the  use  of  the  approximate  energy  ratio  ^  0*8* 

The  conditions  at  the  bubble  minimum,  its  size  and  shape  as  well  as  its 
energy  of  translation  are  not  covered  in  this  analysis*  However,  it  is 
possible  to  calculate  the  position  of  the  bubble  minimum  sinply  by  assuming 
it  halfway  between  successive  bubble  maxima* 

According  to  the  present  state  of  knowledge  the  energy  curves  are  appli** 
cable  to  most  explosives*  Some  deviations  may  be  expected  for  deep  explosions, 
where  only  slight  migration  takes  place*  The  evidence  of  non-migrating  bubbles, 
reference  (j),  shows  a  slight  trsn^  in  the  period  ratios  for  explosives  of 
different  conpositlons*  But,  the  results  of  the  and  C-series  give  no  indica¬ 
tion  of  such  a  tread*  If  desired  the  effect  of  different  explosives  can  be 
accounted  for  in  Figure  6  by  constructing  a  new  curve  which  passes  thrcu£^  the 
appropriate  valtie  at  ^  0  and  subsequently  merges  into  the  old  curve*  But, 

in  most  cases,  such  a  precaution  will  not  be  necessary* 


The  practical  application  of  the  information  obtained  here,  is  discussed 
in  NOI/ER  62-18^4-*  In  this  paper  graphs  are  presented  which  pemlt  a  coavanient 
reading  of  the  bubble  parameters  for  a  widt  ran^e  of  conditions* 
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FIG.  I  EVIDENCE  ON  BUBBLE  MIGRATION 
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FIG. 2  SURFACE  CORRECTION  OF  THE  FIRST  BUBBLE  PERIOD. 

DATA  ARE  FROM  CHARGES  OF  1580  LB  TNT  EQUIVALENT. 
THE  DASHED  LINE  BETWEEN  THE  POINTS  t  INDICATES 
THE  UNCERTAINTY  IN  THE  READING  OF  THE  PERIOD. 
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WATER  SURFACE 


F16.  3a  POSITIONS  AND  RADII  OF  THE  BUBBLE  MAXIMA  CALCULATED  FOR 

THREE  VALUES  OF  THE  MIGRATION  RARAME'^ER  C-  CHARGE  WEIGHT 
1580  LB  -TNT  EQUIVALENT,  FIRING  DEPTH  130  FEET. 
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FIG.  3b  POSITIONS  AND  RADII  OF  THE  BUBBLE  MAXIMA  CALCULATED  FOR 

THREE  VALUES  OF  THE  vilGRATION  PARAMETER  C.  CHARGE  WEIGHT 
1580  LB  -TNT  EQUIVALENT,  FIRING  DEPTH  140  FEET., 
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FIG.  3c  POSITIONS  AND  Rt  j\\  OF  THE  BUBBLE  MAXIMA  CALCULATED  FOR 

THREE  VALUES  OF  THE  MIGRATION  PARAMETER  C.  CHARGE  WEIGHT 
1580  LB -TNT  EQUIVALENT,  FIRING  DEPTH  150  FEET. 
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RATIO  To/T|  AND  D, /Z 
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REDUCED  HYDROSTATIC  HEAD 


FIG  4  PERIOD  RATIO  Jz/lf  AND  RATIO  OF  FIRING  DEPTH  TO  HYDROSTATIC 
HEAD.  ABSCISSA  IS  THE  HYDRO.  TATiC  HEAD  REDUCED  BY  THE 
FIRSI  MAXIMUM  BUBBLE  RADIUS. 
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RATIO  T3/T: 
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FIG.  5  PERIOD  RATIO  Tj/T,  .  TH^  MAGNITUDES  USED  AS  ABSCISSA  REFER 
TO  THE  FIRST  CYCLE  OF  .’ULSATION 
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HWORD  4lB5 
LIST  OF  SYMBOLS 
o  Maximum  Bubble  Radius  (ft) 

«  Minimum  Bubble  Radiiis  (ft) 

=  Reduced  Maximum  Bubble  Radius  (dim  *<^loiiicfis) 

«  See  equation  (19) 

«  Migration  Coefficient  (ft  /lb  ^  ) 

Migration  Ccefficienl;  (dimensionless) 

=  Depth  of  Bubble  Center  at  Bubble  Maximum  (ft) 

»  Firing  Depth  (ft) 

Froude  Number  (Dimensionless) 

=  Surface  Correction  Function  (dimensionless) 

=  Acceleration  of  Gravity  {ft/sec  ) 

=  Distance  to  Bottcan  (ft) 

=  Radius  Constant  (ft^^^/lb^^^) 

=  Period  Constant  (sec  ft^^^/lb^^^) 
a  Subscript  Referring  to  Cycle  of  Pulsation 
a  Chemical  Energy  per  Unit  Mass  of  Eaqolosive  (cal/lb) 
a  Energy  Fraction  Referring  to  Bubble  (dimensionless) 
a  Bubble  Energy  (cal) 

=  Period  (sec) 

=  Reduced  Period  vdim-insionless) 
a  Charge  Wei^t  (lb) 
a  ( D.HO/(  Wi )  (  dimensionless  ) 
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a 

e 


=  Total  Hydrostatic  Head  =  D  +  33  ft  for  sea  water 
»  Total  Hydrostatic  Head  at  firing  depth  (ft) 
m  Surface  Correction  Factor  (dimensionless) 

ts 
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